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Abstract. If the expansion of the early Universe was not purely de Sitter, the statistical imprints of 
■ the primordial density perturbation on the cosmic microwave background can be quite different from 

Q-f 

those following slow-roll inflation. In this paper we study the inflationary signatures of all single-field 
models not plagued by ghost-like instabilities. We assume small deviations from exact scale-invariance, 
as supported by current cosmological constraints, allow for a rapid change of the Hubble parameter and 
the phase speed of scalar fluctuations. We obtain the propagator of scalar fluctuations and compute the 
bispectrum, keeping next-order corrections proportional to the deviation of the spectral index from unity. 
These theories offer an explicit example where the shape and scale dependences of the bispectrum are 
highly non-trivial for reasonable breaking of slow-roll. 
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1 Introduction 



In recent years our understanding of the early Universe has become much clearer as we gain access 
to precise measurements of the cosmic microwave background radiation (CMBR) [1-3]. It is widely 
believed that in the early Universe a period of inflationary expansion took place with scales being pushed 
outside the observable horizon [4]. With this simple idea the Universe becomes very smooth, every 
number density in relics is diluted away [5], and the standard cosmological problems are ameliorated. 
At the same time, distinctive signatures are imprinted in the CMBR anisotropics which can be observed 
by high resolution surveyors, such as the Planck satellite. 

From CMB and large scale structure (LSS) data [6] we know that the spectrum of scalar pertur- 
bations is quasi-scale invariant, with the power being only slightly stronger at larger scales (accord- 
ing to the most recent WMAP7 data analysis, the spectral index of density perturbations evaluated at 
k = 0.002Mpc _1 is n s = 0.969±0.012 [2]). However, a scale invariant spectrum is a model-independent 
feature of inflation. To look for distinctive signatures of different models it is necessary to study higher- 
point statistics of the primordial curvature perturbation, known as non-gaussianities [7], which are a 
result of the inflaton self-interactions and couplings to the gravity sector. 

Amongst the most important of these statistics and potentially easiest to observe is the bispectrum 
or, in quantum field theory language, the three-point correlation function. The first studies on the bis- 
pectrum started with Maldacena's pioneering work [8] devoted to canonical single field-models, later 
generalized by Seery & Lidsey [9] and Chen et al. [10] for theories with a varying speed of sound for 
scalar fluctuations. 1 The latter become known as P{X, cp) models, because the Lagrangian is an arbitrary 
function, P, of the field profile, 0, and its first derivatives through X = — (V0) 2 . These models have 
the property of allowing for large non-gaussianities whenever the speed of sound for perturbations is 
small. On the other hand, models with canonical kinetic terms, even if including higher powers of single- 
derivative operators (suppressed by powers of some high energy physics scale), predict observationally 
small non-gaussianities [12]. For more recent works see, for example, Refs. [13-21]. To understand 
better the microphysics processes which operated in the early Universe one needs to fully understand 
these non-gaussian signatures, since ultimately they have the power to eliminate models against obser- 
vations. With Planck data soon to become available our best hope to achieve this lies within a non-zero 
measurement of the bispectrum, which will be the object of interest in this paper. 

The simplest inflationary models are those where the scalar field theory contains only one scalar 
degree of freedom responsible for sourcing inflation and seeding perturbations: the inflaton field. In 
spite of their simplicity, these models accommodate quite a remarkable variety of interesting scenarios: 
from Dirac-Born-Infeld (DBF) inflation [22, 23] and galileon inflation [17], which are the only known 
radiatively stable theories, to models with canonical kinetic terms. Recently a number of authors have 
also studied models where the galilean symmetry could be relaxed, and focused on the much milder 
requirement that the Lagrangian preserved unitarity resulting in a sensible quantum field theory [24- 
28]. This way the Lagrangian operators give rise to equations of motion for the scalar field which are at 
most second-order in derivatives. Shortly after, these theories were shown to be equivalent to Horndeski 

1 For simplicity we refer to the phase speed of perturbations [11] as the speed of sound. 
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models [26, 28]. First derived over thirty years ago [29], Horndeski theories are described by the most 
generic action of single-field models which do not contain ghost-like instabilities, hence encapsulating 
all the single-field models of interest. 

Until recently the traditional methodology to compute the bispectrum started from the action for 
the comoving curvature perturbation, £ [30, 31], for each model independently. This perturbation is 
the appropriate quantity to study since it is conserved on super-horizon scales, if isocurvature modes are 
absent. However, recent developments have simplified this procedure. First, Gao & Steer [28] and De 
Felice & Tsujikawa [32] (see also Ref. [33]) obtained the universal action for perturbations in stable 
single-field models involving what we call Hornesdeski operators (cubic operators in Q . The only model 
dependent features in the action would reside in different coefficients of each operator. Second, it was 
shown in Ref. [33] that the cubic action for £ has a minimal representation in terms of only five of these 
Horndeski operators, initially derived in Ref. [34] . 

With these recent developments we arrive at a universal methodology to compute the bispectrum 
of single-field models, with the cubic action for £ being always of the form 



;(3)_ ' 



d 3 xdx {aAxC 3 + a 2 A 2 U 2 + a 2 A^{dQ 2 + a 2 h^d&dKd^Q + a 2 A 5 3 2 C(3^- 2 2 } , (1.1) 



where a is the scale factor, dotted quantities are differentiated with respect to conformal (not cosmologi- 
cal) time and d denotes spatial partial derivative. This action appeared in this form in Ref. [34] although 
with different coefficients A ; . 2 As we mentioned before, the model-dependent imprints will be encoded 
in each of the five coefficients A; of the Horndeski operators. There is a priori no hierarchy between 
these coefficients, although specialization to different models can impose specific ratios between A; (DBI 
inflation being one example). The action above will be our starting point in computing the bispectrum. 

Many authors have focused on the study of the action (1.1) under the assumption of slow-roll 
regime in which the inflationary expansion is quasi-de Sitter, whilst allowing for small variation of the 
sound speed of perturbations. In this approximation e = — dlnH/diV, 17 = dlne/dJV and s = dlncJdN 
all obey s, |rj|, \s\ <C 1, where N stands for the number of e-folds. Examples of such works include, for 
example, Refs.[9, 10, 17, 21, 34]. But, what if inflation was not almost de Sitter, and the parameters 
above are not perturbative? Khoury & Piazza [35] showed that this scenario would still produce a 
scale-invariant spectrum, in fair agreement with observations [2], provided the relation 5 = —2s was 
satisfied. 3 We know that both e and 1 17 | must be small to allow for a successful period of inflation. 
However, how small are they required to be? If e is small but not necessarily much smaller than 1, 
then a calculation beyond slow-roll is technically required especially to fit in the current era of precision 
Cosmology. Previous works which have attempted to study correlations beyond the slow-roll regime 
include corrections to the power spectrum in canonical models worked out by Steward & Lyth [39], and 
the work by Gong & Stewart [40] (see also Ref. [41]) who set up a formalism using the Green's function 



2 This action had appeared first in Refs. [9, 10] specialized for P(X, <f>~) models. There the action involved a larger number 
of cubic operators and also a field redefinition. The equivalence of the actions therein and in Eq. (1.1) was explained in Ref. 
[34]. 

3 These authors have also considered solutions within the ekpyrotic mechanism (see Ref. [36] for a review) - see also Refs. 
[37, 38]. In this paper however we focus on the inflationary scenario. 
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method to obtain the propagator of scalar fluctuations to next-next-order in slow-roll. This was later 
generalized by Wei et al. [42] who considered models with a varying speed of sound of perturbations. 
See also Ref. [43] for a calculation report by Bartolo et al. on the power spectrum beyond leading order 
in the context of effective field theories of inflation [44, 45] . 

In this paper we will be interested in models where the spectrum is almost scale-invariant, so that 
the condition found by Khoury & Piazza is mildly broken and becomes s = — 2e + 5. This requires working 
to all orders in e and 5, but perturbatively in 5, in deviations of the spectral index from unity, n s -Kl, 
and time variations of e and s. This last assumption ensures that the conditions s, \s\ < 1 are preserved for 
a sufficiently large number of e-folds, which is a primary requisite for inflation to last at least 60 e-folds. 4 
In this sense, this is a calculation beyond the (strong) slow-roll approximation. We refer to it as a scale 
invariant approximation. Starting from the action (1.1) the calculation is immediately applicable to all 
Horndeski models. For this reason we will work with arbitrary interaction coefficients A; - assignments 
to these coefficients will correspond to specific models. We organise the calculation in increasing powers 
in the hierarchy of slow-variation parameters, focusing on leading and next-order contributions. This 
organisation scheme follows the one used throughout Refs. [10, 17, 34] and generalizes their results. 
We explicitly obtain the scale dependence of the bispectrum, which appears more intricate than the 
one found with the assumptions of (strong) slow-roll. In particular, the bispectrum exhibits a strong, 
power-law scale dependence, accompanied by the traditional weak, logarithmic dependence. 

Outline. — This paper is organised as follows. In §2 we obtain the dynamical behaviour of the scale 
factor and the other relevant background quantities. The dynamical behaviour is exact in e and s, and 
perturbative in the time variation of these parameters and n s — 1. In §2.1 we obtain the power spectrum 
for scalar fluctuations in the quasi scale invariant approximation, and we derive in §2.2 formulas for the 
elementary wavefunctions from which the scalar propagator is built to next-order in scale-invariance. We 
use the results obtained in the previous sections to calculate the bispectrum of perturbations in §3 and 
comment on the differences when using the slow-roll approximation. We conclude in §4. The appendices 
collect useful formulas and detailed expressions to recurring integrals necessary to produced closed-form 
bispectra. 

Notation. — We choose natural units in which c — h — 1 and we take the reduced Planck mass, M P = 
(87iG) -1 / 2 , to be unity. The metric signature is (—,+,+,+), and purely spatial indices are denoted by 
Latin letters {i,j, ...}, which are also occasionally used to refer to each Horndeski operator. Derivatives 
with respect to conformal time, t = dt/a(t), will be denoted by dots. 

2 Background evolution beyond exact scale invariance 

We are interested in quasi-scale invariant models which fall under the class of Horndeski theories. As 
argued in Ref. [35] a scale invariant spectrum for perturbations can be accommodated by a background 

4 This work explicitly excludes models with features in the potentials which can trigger the slow-roll parameters to temporar- 
ily grow during inflation (see for example the review by Chen [46]. 
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where both the Hubble parameter and the speed of sound of perturbations vary significantly in time 
provided they obey the exact relation s = —2s. 

We start by considering an arbitrary Horndeski theory for a homogeneous scalar field, on which 
small, inhomogeneous perturbations, 5cj), develop. All we require is the spectrum of perturbations 
to be quasi-scale invariant. The comoving curvature perturbation is related to these perturbations via 
£ = aH5cj)/(j). Often, one computes the correlations functions in the comoving or £-gauge rather than 
in the uniform-density or 50-gauge. This has an important advantage: the correlation functions of £ 
will be time-independent whereas the ones for 5cj) will not. For this reason it is more convenient to work 
with the action (1.1) since the primordial comoving perturbation, £. We will come back to this point 
shortly. 

For all inflation models involving one single clock in the Universe, the quadratic action for the 
primordial comoving perturbation, £, can be written as 



S™ = - 
2 



d xdza % 



(2.1) 



where z is required to be a well defined, differentiable function, but otherwise arbitrary. 5 The time 
evolution of z(y) will be parametrized by w = dlnz/diV. For consistency, we also work to all orders 
in w. For an arbitrary w, Khoury & Piazza's relation between e and s required for scale invariance, is 
generalized to 6 

3s = — 2e — w . 

This formula reduces to s = — 2e for P(X,4>) models. It is more convenient to study the dynamics of 
background quantities as a function of y, satisfying dy = c s dr, rather than conformal time, t. 7 In this 
new time-coordinate, the action (2.1) then becomes 



:(2) _ 



d 3 xdyq 2 C /2 -(9C) : 



(2.2) 



with q = dy/zc s ' and where the prime denotes derivative with respect to y. Introducing the canonically 
normalised field v k = q^ k , where the subscript k indicates the Fourier mode, and changing variables to 
v k = A fc y—ky, we find that the modes A k obey a Bessel equation of the form 



e-{ q - + 



4y' 



= o 



The solution to this equation is a linear combination of Hankel functions, H^\—ky) and H^ 2 \—ky), 
with order 



3 

v — — 
2 



n c — 1 



5 For P(X,cp) models, z = s/cf [34], but the same need not be true for other models. That is the case for galileon inflation 
theories Ref. [17], for example. We will come back to this point in §2.2. We note that z here is not the usual z for fc-inflation. 



6 We refer the reader to Eq. (A.3) in appendix A for details. 
7 The reasons for this choice will be clarified in appendix B. 
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For perfectly scale invariant spectrum of perturbations then v = 3/2 and the spectral index is unity. This 
results in the power being constant in all scales. See appendix A for details of how the spectral index 
relates to the background quantities - the precise formula will be unimportant in what follows. 

Agreement with the appropriate normalization of the propagator gives the overall evolution of the 
elementary wavefunction as follows 



IP J-ky 



(2.3) 



8k a^Jzc s 

Finally, the power spectrum of perturbations can be obtained by evaluating the scalar propagator at equal 
times, yielding 

C-y) 



(ak 1 )«k 2 )) = (27r) 3 5^(k 1 + k 2 ) 



71 



8 a 2 iy)z(y)c s (y) 



H^X-ky) 



(2.4) 



To further simplify this expression we will require the evolution of the scale factor a in the y coordinates. 
In order to do so, we specify the parameters in the slow-variation catalogue, which will allow us to 
perform a uniform expansion up to next-order terms in the quasi scale invariant approximation. 



Slow -variation catalogue. — To obtain the dynamical evolution of the scale factor in a background with e 
and s smaller than 1 (but not necessarily much smaller than 1), we introduce the following slow- variation 
parameters: 

dine dins 
n = — — and t = — — . (2.5) 
' dN dJV 

We assume these parameters satisfy \rj\, \t\ <C 1. 

To determine the behaviour of the scale factor, we start by integrating c s d.T to compute y. Working 
perturbatively in r\ and t but to all orders in e and s, we find 

+oo m 

J = — 77 J] {(* + s) m + (e + s)™- 1 £ fc(ts + erj) } . 

m=0 k=0 

These converge and give the dynamical behaviour of a in y-time as follows 



a(y) = - 



1 + 



sr] + ts 



(1 



■s) 2 J 



(2.6) 



Hy(l-e-s) L 

This result was first obtained by Khoury & Piazza in Ref. [35]. Using Eq.(2.6) we can further simplify 
the two-point correlator (2.4) in y-space, which becomes 

7iH\y){-yf{l-e-sf 



(C(k 1 )ak 2 )) = (27i) 3 5^(k 1 +k 2 ) 



8 



z(y)cKy) 



1-2- 



£17 + ts 



'(1 



(2.7) 



In a generic model, the two-point correlator will evolve in time. However, for a single-field model, 
by virtue of conservation of the primordial curvature perturbation on super-horizon scales [47, 48], 
it follows that the power spectrum will rapidly converge to its asymptotic value. 8 Therefore, we can 



8 See Ref. [49] for comments on the number of e-folds necessary for such asymptotic behaviour to be reached after horizon 
crossing. There the authors study a two-field inflation model, but the same conclusions still apply in single-field models. 
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focus our calculation a few e-folds after horizon crossing for a given scale, k+. This has become the 
standard approach when calculating the asymptotics of the power spectrum. This remains true even in 
multi-field inflation, where in this case this procedure sets the appropriate initial conditions for evolving 
correlation functions after horizon crossing, using techniques such as 5N [50], transport equations [51, 
52] or transfer matrices [53] (see, for example, [54—59]). The advantages of applying this approach 
are twofold: on one hand, at this point in the evolution the elementary wavefunction is completely 
characterized by the growing mode, as the decaying mode has become negligible; on the other hand, 
the power spectrum obtained in this way is already the asymptotic value, which corresponds to the 
late times regime. In other words, at this point of evaluation the primordial curvature perturbation has 
already classicalized. 

With this methodology in mind, and working perturbatively in 17 and t we find that 

17* 



1- 



SKS+ 



1-E+-S 



■H-Ky) 
W-Ky) 



(2.8a) 
(2.8b) 



to next-order in our approximation. In the remainder of this paper we employ the same notation as 
in Refs. [17, 34] where starred quantities are to be evaluated when some scale, k+, has exited the y- 
horizon, that is, when k+y+ — — 1. It is important to leave this scale arbitrary since this will produce 
the scale dependence of the correlation functions. The appearance of the logarithmic term corresponds 
precisely to the number of e-folds, N+, in the y-coordinates, which have elapsed since the scale k+ has 
exited the horizon. 9 The expansion in Eqs. (2.8) is only valid up to a few e-folds outside the horizon, 
when one can trust the Taylor series expansion to first order. We note that quasi scale invariance demands 
17* = t*, up to next-next-order corrections, which are beyond the scope of this paper. 

Given the individual dynamics of c s (y) and H(y) deduced in Eqs. (C.3) and (C.5) of the appendix 
C, we can replace for the evolution of the scale factor in Eq. (2.6). We find: 



a(y) = 



CsiX-Ky) 1 - < *-* 

H*(l-e*-sJ(-y) 



H-Ky)+ 



2(l-e*-sJ 



(ln(-^y)) 2 , (2.9) 



where the next-order parameters a and /3 are defined in Eqs. (C.2). Despite its complicated structure, 
it correctly reproduces the results of Refs. [34, 35] when we take the limit of exact scale invariance, 
for which 17* = = 0. The extra terms in Eq. (2.9) are precisely the corrections to a pure power-law 
evolution, since they correspond to the perturbative expansion of the dominant power law behaviour. 

2.1 The scalar power spectrum 

As mentioned above, the power spectrum, defined as the propagator evaluated a common time 

(OkOftka)) = (2ti) 3 5 (3) (k! + k 2 )P(fc) (2.10) 



9 The notation for the number of e-folds as a measure of expansion was initially defined by Sasaki and Tanaka in Ref. [60] . 
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where k = |k| = IkJ = |k 2 |, will in general depend on time. Taking the super-horizon limit, \ky\ — * 0, we 
find that the dominant contribution arising from the Hankel function in Eq. (2.7) is given by 

H^i-ky) — --p-j {l + („ s - 1) [ - 2 + rE + ln(-2fey)] } , 

where y E ~ 0.577 is the Euler-Mascheroni constant. Whatever y -evolution the remaining terms in Eq. 
(2.7) have in this limit, they should precisely cancel the y dependence of the Hankel function. This is a 
requirement imposed by conservation of £ on super-horizon scales in single-field models. 10 This allows 
us to write the power spectrum as 

P(k)= \ * 3 * {l + 2E* + (n 5 -l)ln(fc/fcj} , (2.11) 

where we have defined 

E* = -&+^— [-2 + r E + ln2] . (2.12) 

This is the power spectrum for a quasi-scale invariant Horndeski theory. To recover the perfectly scale 
invariant formula obtained by Khoury & Piazza [35] we simply set and n s — 1 to vanish, for these enter 
the power spectrum as next-order corrections only. Whilst being time-independent, Eq. (2.11) correctly 
reproduces the expected logarithmic scale-dependence in the strong slow-roll regime, in which e and s 
are taken to be slow-roll parameters. Indeed, we can check that defining the dimensionless version of 
the power spectrum using the usual rule & = k 3 P(k)/2n 2 , this formula agrees with 

dln^> 
*- 1 = dTnT ■ 

Again, the explicit expression for n s — 1 will not be required in any stage of the calculation. It suffices to 
treat it as an arbitrary next-order quantity. 

2.2 Obtaining next-order corrections to the wavefunctions, £ k 

We are interested now in obtaining the bispectrum of perturbations starting from the action Eq.(l.l). We 
will apply the same procedure defined in the previous sections and rely on the slow-variation catalogue 
introduced there. Since we will be keeping up to next-order corrections in the slow-variation parameters, 
we also require corrections to the elementary wavefunctions up to the same order. It is easier to think 
about these corrections by locating them in a Feynman diagram corresponding to the three-point function 
in question. The sources of next-order corrections are therefore quite well-defined and we explore them 
in what follows. Next-order corrections were first enumerated by Chen et al. [10] for P{X, 4>) models. 

External legs. — The corrections to the external lines of Feynman diagram correspond to evaluating the 
propagator in the asymptotic regime, when \ky\ — * 0. From Eq.(2.11) we can read off the external legs 
corrections and conclude that the propagator becomes 

r „ n i HJ l — e^ — s^r — 1 i 

C fc (y) (eJCter " a ° A J 1 + Ht + _L- h(t/t . ) . (2 .i3) 

2 y /zT(kc sic y> z <- 2 J 



10 The time-independence of the power spectrum in single-field inflation was carefully explained in Ref. [34] to which we 
refer the reader for more details. 
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As argued in §2.1 the time independence of this value is guaranteed because the primordial perturbation 
£ is conserved on super-horizon scales. 

From the requirement of time independence of the power spectrum it also follows that whatever 
function we assign to z(y) it must generically behave, at leading order, as 



~ (.-Ky) 1 -'*- 



(2.14) 



This represents a mild requirement on an otherwise arbitrary but smooth function z(y). We check that 
for P(X, <fi) models, where the formula for z is well-defined, z = e/c 2 , this behaviour is obeyed by virtue 
of Eqs. (2.8) and (C.3), together with the condition s + = — 2e* (which holds at leading order in our 
approximation scheme). 



Internal legs. — This correction is a result of the spectrum being slightly tilted. This manifests directly in 
Eq. (2.3) for the elementary wavefunction which, using Eq. (2.6), we rewrite for clarity as 



,V m>)(-jWi-<-i) l, sn + a x m 



(2.15) 



Since the order of the Hankel function ,v, differs from 3/2 by next-order corrections, these originate in 
turn corrections to £ fc . To evaluate these we start by Taylor expanding the Hankel function around order 
v = 3/2. Following Refs. [10, 17], we find that the background evolution of the wavefunctions is given 
by 

„ , „ i// (1 - e -s ) 

(2.16) 



Cfc(y) (bacfc§ro "" d) = — £ * *?h l - iky)e iky 



which agrees with the leading order behaviour of Eq. (2.13) in the limit \ky\ — » 0. The next-order 
corrections arising from 5v = v — 3/2 = — (n s — 1)/2 and the slow-variation of the remaining background 
quantities in (2.15) organise themselves in what we label as internal leg corrections as follows: 

a t M \H. (1 — £± — S+) 
^internal) * v * * J 



2yzT(fcc s J 3 / 2 



e iky {l-iky) 



+ 



n c — 1 



e~ iky {l + iky) 



-dC-2e 



iky 



TZ 



i_ e i*y (1 _ ifcy) 



(2.17) 



The integral representation written above corresponds to the exponential integral function, Ei(— ky), 
defined for real, non-zero argument: 

moo 



rx e* 
— dt 
t 

00 



The time dependence of the internal leg corrections is fairly intricate and verifying that there are no 
divergences in Eq. (2.17) in the late time limit represents a minimal check on this result. We indeed 
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confirm this is the case. Expressions for time derivatives of the wavefunctions can be found in appendix 
C. 



Corrections arising from interaction vertices. — The interaction coefficients A ; can generically evolve very 
fast, and one expects 

A i ^A i ,(-/cy)"|l + aln(-/c,y) + b(ln(-/c,y)) 2 J , (2.18) 

with the power n depending on the exact expression of A iis and a and b being next-order terms. This 
power law dependence needs to be taken into account for consistency with the scale invariant approx- 
imation, and we will study an explicit example in §3.2. There we explain in detail that this power law 
behaviour can potentially bring problems to the convergence of the integral. In particular, we might have 
to require that only some values of e are allowed in order that the three-point correlator does not evolve 
in time. Mainly to simplify our results, we assume the interaction vertices to have a very smooth and 
slow evolution in the y-coordinates. To this end we introduce a supplementary slow-variation parameter 

1 dA ; 

hi = — -— - (2.19) 
AjH dt 

which satisfies <C 1. We expect this approximation to be valid in all Horndeski models where the 
interaction vertices are slowly varying in time, even though the respective assignments for A ; can be 
rather intricate. As described in §3.2 the methods of this paper can still be used to compute the n-point 
correlation functions in models where this assumption fails. We conclude that, if A ; is slowly varying, 
then each interaction vertex evolves as 



A;~A 



i-- — - — H-hy) 

l-e^-s* 



(2.20) 



This means that for simplicity in the remainder of the paper we will work with n = in Eq. (2.18), but 
we will show how to deal with n ^ in one simple example. This assumption concludes the presentation 
of the slow- variation catalogue. 

We note that a further simplification to the background dynamics follows from inspection of Eqs. 
(2.13) and (2.17). The variable w that parametrizes the variation of z, does not appear in our formulas. 
Indeed such term is only present implicitly since it is encapsulated in the spectral index n s — 1 . Therefore 
the calculation is independent of our choice of parameter w. As a consequence we can effectively reduce 
the relation we found between e, s and w, to one which only involves the first two variables. To this end, 
we take the relation between e and s to be 

5 = -2e + 5 

where 5 is a next-order parameter. We emphasize that this procedure in no way lacks generality: the 
calculation is still non-perturbative in e and s. The formula above generalizes the regime studied in 
Ref. [35] by considering a quasi-scale invariant spectrum. This also implies that t = rj, which is in 
agreement with the observations in §2. This procedure dramatically simplifies our formulas and reduces 
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the slow- variation catalogue to a minimum of parameters, which we choose to be 5 and 17. The quantum 
field theory calculation we will present is thus only non-perturbative in the e parameter. 11 We note that 
this will imply some conceptual restructuring of our formulas. In particular, Eq. (2.16) now contains 
background and next-order contributions in 5+, which can now be interpreted as an extra contribution 
to the external legs at next-order in Eq. (2.13). 

Eqs. (2.13), (2.17), and (2.20) combine to produce the overall corrections which are required to 
produce the bispectrum at next-order. They naturally combine to produce two separate contributions 
to the bispectrum: the background contributions and next-order corrections arising from the vertex and 
external legs; and the next-order contributions arising from the internal legs of the propagator. We will 
use this way of partitioning the bispectrum when presenting the results for the bispectrum in §3.1 by 
labelling the first contributions as type a bispectrum and type b bispectrum for the second. 

3 Non-gaussianity in the bispectrum 

Non-gaussian signatures encoded in the CMBR work as a powerful discriminant between inflationary 
models. Quantum field theory correlation functions were initially studied by Schwinger [61] and Keldysh 
[62] in the 60s, and later applied to Cosmology by Jordan [63] and Calzetta & Hu [64] . But it was only 
with Maldacena's publication in 2002 [8] that the applications of the in-in formalism to non-gaussianity 
were made more clear, together with a pair of papers by Weinberg [65, 66]). This formalism is the 
appropriate construction to compute expectation values, and we briefly summarize its basic ideas (there 
are a number of papers which review the in-in formalism - see, for example, Refs. [67, 68]). 

In this paper, we will be interested in computing three-point correlation functions, B, at tree-level in 
the interactions of the comoving curvature perturbation, £. In terms of the Argand diagram in complex 
time y, these correlations are obtained by performing a path integral from the true vacuum of the theory, 
at t — > — 00, to the time of interest when we compute the expectation value, added by the path integral 
performed backwards to —00. Schematically, this can be translated into 

(C(k x )C(k 2 )C(k 3 )) = J [dC + dC_] C + (ka)C + (k 2 )C + (k 3 ) e is ^~ is ^ (3.1) 

where the forward path integral is labelled by the fields £ + whereas the backwards path integral is 
labelled by the fields £_. 12 These fields are constrained to agree at any one time later than that of the 
observation, such that the path integrals translate into a contour of integration turning and crossing the 
real y axis. From this construction we see that we need the cubic action, S, for the perturbation £, 
written in Eq. (1.1), to obtain the lowest order non-vanishing Wick contractions between the fields. This 
is because at leading order in the interaction, the first non-vanishing contribution arises from contracting 

11 Our calculation is non-perturbative in e and s. By realizing that we can write s = —2s + 5, this allows us to replace the 
non-perturbative parameter s by the perturbative parameter 5. In the remainder of the text, the dependence in 5 is absorbed 
into 5. 

12 In practice, to project onto the true vacuum of the interacting theory, one needs to translate the integration contour to 
slightly above and below the negative real axis of conformal time, respectively. This prescription is in many ways similar to the 
is trick recurrent in quantum field theory and ensures the convergence of the integral. 
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the three fields inside the correlator (essentially three copies of £) with the three fields in the cubic action. 
As the action is composed by five Horndeski operators, the total bispectrum will be given by the sum of 
the corresponding five (individual) bispectra. The form of Eq. (3.1) also ensures that the contribution 
from the backwards path integral is precisely the complex conjugate of that of the forward integral. 13 
Therefore, it suffices to compute one of these path integrals. 

There are several ways of presenting the bispectra. Some authors prefer to present their answers in 
terms of a rescaled bispectrum traditionally defined as [50, 69] 

6 P{k l )P{k 2 ) + cyclic permutations 

Cosmological constraints for / NL are quoted in Refs. [70, 71] . When calculating / NL we need to add two 
symmetrized copies of the power spectrum - this is what we mean by "cyclic permutations". To place 
constraints in model building, it is common to write / NL evaluated in some specific momentum configu- 
ration, such as the equilateral limit. We then compare this estimate with the amplitude of the bispectrum 
in the equilateral template used in CMB analysis. However, simply because some bispectrum shape peaks 
in the equilateral limit, does not mean it perfectly matches the bispectrum equilateral template. Hence, 
in order to place constraints on the parameters of some theory we ought to further require the evaluation 
of the correlation between the bispectrum we are studying and its closest template, since it is desirable 
that the error bar associated with the constraints is comparable to the resolution of the data. 

Others have observed that the entire bispectrum shape can be used as the individual signature of 
each inflationary model [72, 73]. Decoding the entire bispectrum shape into primitive or fundamental 
harmonic shapes by performing a partial-wave expansion acts in the same way as the decomposition 
of the angular power spectrum in terms of Q, once we identify the number of the harmonic with the 
multipole moment, I. Whichever final bispectrum shape is produced it will be a linear superposition 
of five individual bispectra, with weights given by the respective interaction coefficients, Aj. This has 
a major significance for the overall bispectrum shape: it can be viewed as a unique fingerprint of the 
Lagrangian structure of the theory. In this perspective, there is no real advantage in specifying the results 
for / NL , when one can use the bispectrum shape as the most sensitive and discriminating degree of 
freedom. Which of these methods will prove to be more efficient in data analysis is still not absolutely 
clear, although the modal decomposition seems to be quite promising (see, for example, Ref. [74] for 
comments on the efficiency of this algorithm). In particular, it is possible to derive consistency relations 
between the amplitudes of each harmonic, which make it easy to rule out classes of inflation models [21] 
whenever these relations are not supported by data. 

In this paper we shall adopt this last philosophy and present the general formulas for the bispec- 
trum. The study of shapes is however beyond the scope of this paper, and we leave it to forthcoming 
work. We specify estimates of / NL in specific applications only. 



3 This is only true at tree-level for diagrams which involve one interaction vertex only. 



-11- 



3.1 Bispectrum of Horndeski theories 

As argued before, to compute the non-gaussianities it will be convenient to write the action for pertur- 
bations in y coordinates. From Eq. (1.1) this is given by 



:(3) 



c 2 1 

d 3 xdya 2 \^A 1 C 3 + c s A 2 ^ ,2 + -A3ad0 2 + c s A 4 Cd i i:d l (d- 2 n + c s A 5 d 2 ad i d- 2 a 2 
a c s 

(3.3) 



Our calculation is organised as follows: 

Bispectrum type a. — The background contributions to the bispectrum and the corrections arising from 
the slow-variation of the interaction vertices and external legs can be written in the form 

B a = Al * 6 H * ( At E *l N^k,) { Re (P a (fc x )j r ) + Re (o a {k 1 )J r+1 ') + T a (kA + cyclic permutations , 

(3.4) 

where the addition of cyclic permutations entails exchange of momenta fc x — » k 2 — * k 3 . The functions 
N a (k 1 \ P a (k 1 \ Q a (fc 1 ) and T a (k 1 ) are listed on table 1 and table 2 contains the assignments for the 
functions jL 

We note that the coefficients on table 2 for the operators ^' 2 , gd^dfi' 2 ?, and d 2 ^^' 2 ?) 2 
are equivalent. Having done the calculation independently for each operator, this is a minimal check of 
the correctness of our results. This is because these operators have the same time derivative structure, 
and differ only in the arrangement of spatial derivatives, and therefore, in the momentum dependence. 

Bispectrum type b. — Likewise, the contributions from the propagator in the internal lines of the Feynman 
diagram can be consolidated for all operators in the form 

Bb = ^wn.fcj **fco{(n, - d£ [ Re (p/^w+Re (QKwy^ikS 

+ Re (i? b (fc!) J r ) + Re (s 6 (fc x ) J J+1 ) + (n s - l)T b (fc x ) j (3 ' 5) 
+ cyclic permutations , 

where the addition of cyclic permutations entails exchange of momenta fc x — * k 2 — * k 3 . The functions 
N b (k 1 ), P ; & (fci), Q 6 (fci), i? b (fci), S b (k 1 ) and T 6 ^) are listed on tables 3 and 4. Table 5 contains the 
assignments for the functions I r and J r 



Both types of corrections are computed using the same in — in formalism rules. The third Horndeski 
operator, £(<5£) 2 , however, presents an additional degree of complexity (which also had to be dealt with 
in Ref. [34]). We briefly review it here. We start by noting that this operator is undifferentiated (with 
respect to time), and therefore the integrand will be, at least, power-law divergent. This is because 
there will be insufficient powers of y, to counteract the powers in denominator owing to the presence 
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of a(y) in the integrand. As a result the integrand behaves, at leading order, as a 2 /c s ~ 1/y 2 . 14 The 
next-order corrections will contribute with logarithmic terms. This has serious repercussions in the final 
result because it can potentially lead to power law and logarithmic divergences in y when one takes the 
limit \y\ — » - these are dangerous interactions, as named by Weinberg [65]. In order to ensure that the 
real part of the correlation function converges in the the asymptotic limit, one needs to carefully isolate 
the primitively divergent contributions: either power laws like |y|~ a , with a > 0, or logarithmic terms, 
ln(— y ). This is done by integrating all divergent integrals by parts. Only two situations occur as a result: 

i. the isolated divergent terms contribute with a purely imaginary part to the correlation function. In 
this case, this divergence becomes irrelevant since the total correlation function is given by the sum 
of two path integrals which are complex conjugates. 

ii. the isolated divergent contributions are real and contribute to the final answer. However when we 
sum type a and type b bispectra we verify that these divergent contributions cancel out amongst 
themselves. The final result contains only finite contributions. 

It is crucial to take into account these two possibilities to obtain a correct result: first, to ensure that 
the correlation functions do not diverge in their asymptotic limit; and second, to guarantee that the 
calculation contains all the relevant convergent contributions to the overall bispectrum. This is particu- 
larly important when checking whether our results are consistent with Maldacena's factorization theorem 
[8, 75] - this asserts that the bispectrum, in the limit when one of the momenta is small, should factor- 
ize into two copies of the power spectrum multiplied by n s — 1; in other words, in this limit, two and 
three-point correlators talk to each other. 

3.2 Features of the bispectrum of Horndeski theories 

These results extend those obtained following a strong slow-roll inflationary phase in the early Universe. 
Some observations following our formulas share the same ideas as in previous works in the literature, 
and we summarize them here for completeness: 

Enhancement of non-gaussinities. — It is apparent from our formulas that the enhancement of the bispec- 
trum can be very different from the one found in models under the slow-roll approximation. This strictly 
depends on the coefficients A ; . Just like observed in Ref. [17], if the interaction vertices in the action 
(1.1) do not contain enough powers in the speed of sound c s to counteract those in the denominator of 
the bispectrum, then the overall dependence can be stronger than that commonly associated with DBI 
models / NL ~ c~ 2 , corresponding to a bispectrum behaving as B ~ c~ 4 . Depending on z, more exotic 
models could reproduce such behaviour. To gauge whether this scenario would be permissible could in- 
volve applying the partial wave decomposition method to the Horndeski overall bispectrum shape. Five 
measurements would be required to fix each of the A ; interactions, and a number of additional mea- 
surements to break the degeneracy in other parameters, such as e and c s , in A ; . Ultimately this would 

14 The calculation of the bispectrum of this operator is very similar to the one presented in Ref. [34], which was restricted 
to the slow-roll approximation. The reason this happens is because, contrary to the remaining operators, the constant y in the 
integrals in Eqs. (D.l) and (D.7) is an integer, which allows us to write compact results in tables 1 and 3. 
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allow us to place constraints on all the parameters of the theory. Only then would we be able to conclude 
whether there would be enhancement of non-gaussianities in these models. 



Presence of logarithms. — Our formulas contain logarithms of momenta, encoded, for example, in the 
master integral J y (defined in Eq. (D.7)). There are two varieties of logarithms as noted in Ref. [34]: 
those which depend on the reference scale, ln(k+/k t ) and those which depend on the perimeter momen- 
tum scale, lnf^i/fc,.). The first type of logarithms are clearly of the same nature as the ones identified in 
the power spectrum (2.11): they encode the scale dependence of the quantities they multiply. One can 
choose the reference scale, k+, to minimize these logs; alternatively, one can use this degree of freedom 
to measure primordial non-gaussianities on different scales. The last type of logarithms are shape depen- 
dent. To better understand this, one can write one of the three momenta, k h in terms of the perimeter, k t , 
and two additional coordinates, describing the angular dependence. This implies that ln(fc ; //c t ) is effec- 
tively only a function of the angular coordinates, and it is therefore responsible for the shape dependence 
to the bispectrum. 

Away from (strong) slow-roll. — If the inflationary background is almost de Sitter, it is easy to see that 

the power law behaviour k+ e * , where a is some integer, can be written as a next-order logarithmic 
correction by performing a Taylor series expansion. This can be checked by taking the limit of very small 
e and comparing our results to those of Ref. [34]. Away from the strong slow-regime, the dependence on 
the reference scale, k+, arises from non-expandable power laws in addition to logarithmic contributions. 
For comparison, let us write the formula of the bispectrum type a for one operator, say £ /3 . We find that: 

5s* 



K ( K 
A^ + ln— B„ + C*ln — 



sin(?)r(?)| 



+ xp {0 \3) ffi, + 2Q In ^ + C^ (0) (£) ) 
fc £ 1^ 



+ cyclic permutations , 

(3.6) 

with E, = 3 + rj^jr and q = ^ ■ This expression exhibits an unusual power law dependence on 
the reference scale, fc*, which is absent from previous studies, where only the logarithmic contributions 
ln(k + /fc t ) were known [34]. By keeping the reference scale k+ arbitrary in our calculation, the scale 
dependence of the bispectrum can be calculated. This will include the contribution of the power law 

_3_ 5g * 

scaling as k t 1+E * , which could receive large corrections from e+. In Ref. [35] Khoury & Piazza chose 
k ir = k t , which masks the power law effect. 

Considering a slow-roll expansion, by which we take to be a slow-roll parameter (treating it on 
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equal footing as 5 and n s — 1), Eq. (3.6) simplifies to 

f<n 3A U H* fh u 3(n.-l) <5„ 

B ?« - tf^tf i t - !+if - 1 + -3— ( " 2 + rE) + t ( - 13 + 6 ^ + t (25 - 

3(n s - 1) 



+ 



f2k 1 k 2 k,\ kA 
ln l^f^J~ (hl * + 35 *~ 5Cjln ^J 



which resembles the bispectrum obtained in Ref. [34]. 15 We note that in Eq. (3.7) the dependence on 
k+ appears through logarithms, and the power law behaviour has explicitly disappeared from the result. 
This is because by taking e to be a small parameter, the power law is well described, at next-order, by 
a logarithmic contribution. We also observe that the reference scale, k+, appears in the form of two 
logarithmic functions: Inikjk^) and ln(k+/k t ). These species of logarithms were thoroughly studied in 
Ref. [34] to which we refer the reader for details. 

The behaviour described above shows that whereas the power spectrum has a universal weak log- 
arithmic scale-dependence, cf. Eq. (2.11), the bispectrum of these theories beyond the slow-roll regime 
can have a much stronger scale-dependence. In principle, this could be used to distinguish, on observa- 
tional grounds, between these theories: from CMB (k < 0.5frMpc _1 ) to cluster scales (k > 0.5frMpc _1 ), 
interpolating with scales probed by the galaxy bispectrum. A number of authors have studied constraints 
on non-gaussianity arising from galaxy surveys [76, 77], including its relation with biasing [78]. The 
scale dependence of the bispectrum was initially studied by Chen in an infrared model of DBI inflation 
[79], and later investigated in P(X,4>) models (a subset of the Horndeski class) by LoVerde et al. [80]. 
It is not the aim of this paper to present a detailed analysis of the scale dependence of the bispectrum. 
We rather want to offer an explicit example of theories which inherit a highly intricate scale dependence 
from the background dynamics. 

Comparison with previous results 

In certain limits, some of our results overlap with others in the literature. Here we compare them. 

Khoury & Piazza. — First, we recall that Khoury & Piazza [35] have calculated the bispectrum of La- 
grangians involving only the inflaton and its first derivative in an exactly scale-invariant background, in 
which s+ = — 2e+. In this paper we have extended this study in two ways: by performing the calculation 
perturbatively away from the exact scale-invariance regime and applying it to the Horndeski class of 
models. 

In particular let us take a specific example for comparison. Focusing on the operator ££ /2 and in 
the action studied by Khoury & Piazza we take 

A 2 (y) = 4( e - 3 + 3c ') • 



15 To be more precise, we can indeed show that our results for the overall bispectrum agree (summing type a and type b), 
including the explicit dependence on the reference scale, fc + . 
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We note that A 2 (y) is rapidly varying in time and should therefore be placed inside the integral to 
compute the bispectrum: 



B^I2 ~ 



2 



dyA 2 (y)c s (y)a z (y)y z {l-iky) , 



where the factor y 2 comes from the two time derivatives of the wavefunctions and the last one from 
the undifferentiated wavefunction. For comparison purposes, we only retain the contributions at leading 
order in scale-invariance, which means setting all the corrections we have focused on this paper to zero. 
Therefore, selecting the leading order contribution to our bispectrum in Eq. (3.4), we find 

+ (e * " 3) (I) ^ r i 1 ~ irt) cos (r^t ) ft + *. + - 3 *.» } 

+ cyclic permutations , 

(3.8) 

where we have retained the scale dependence through k+. This formula reproduces the results of Ref. 
[35] provided we choose the reference scale to satisfy k+ = k t . 

Whenever comparison was possible our results reproduce those of Ref. [35] . 

Noller & Magueijo. — Second, in Ref. [81] Noller & Magueijo estimated the magnitude of the next-order 
corrections to the bispectrum of models whose Lagrangian was a polynomial function of the field profile 
and its first derivatives - the so-called P(X,4>) models. In comparison, our calculation extends to the 
larger Horndeski scalar field theories. Their estimate focused on an approximation of a subset of next- 
order corrections important when |fcy| — > 0. However, there are other contributions which contribute 
equally to the bispectrum and are sensitive to the dynamics around the horizon crossing, when the 
approximation \ky\ — > fails. As we have shown these corrections can be obtained by evaluating the 
change in the Hankel function when we take the order to be v = 3/2 — (n s — l)/2, rather than at the 
exact scale invariant choice, v = 3/2. 

Moreover, in our calculation we have performed a uniform expansion in next-order corrections 
in the quasi-scale invariance approximation, using in particular Eq. (2.9) for evolution of the scale 
factor, a, and Eq. (C.3) for the evolution of the sound speed of perturbations, c s . Ref. [81] assumed 
a perfectly scale-invariant background, on top of which perturbations would develop. This amounts to 
setting a+ = /S v = in our Eq. (2.9), even when their contribution is as relevant as terms linear in n s — 1, 
as we have previously argued. As we have seen, the scalar fluctuations are sensitive to the background 
dynamics, and in particular to the next-order corrections we have calculated. 

For completeness, we will investigate one particular limit studied in Ref. [81] when the speed of 
sound of fluctuations is small, c s <C 1, in P(X,cp) theories (these include DBI inflation). This regime 
is phenomenologically interesting since it is known to be related with large non-gaussianities [23]. To 
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make this application more explicit we apply the notation of Ref. [34] (see also Ref. [10]), and we 
observe that the dominant operator in this regime is £ /3 



;(3) 



d 3 xdyaA lCs 2 C /3 



with 



-c 2 -2c 2 - 

Hc s 4 V s s S 



where 

A lf2f x + l 



£ 6 \ c 



fx is taken to be a constant. We can see that taking the limit of small c s , corresponds to considering 
|-| » 1, with which the interaction vertex is well approximated by 

Inspection of this formula reveals that this interaction vertex is rapidly varying, by virtue of its depen- 
dence in e, H and c s , of which only e is slowly varying. 

Our formulas are easily adapted to the case of rapidly varying A,. For A^y) we find 



Ai(y) = A u (-kjO 



1 + 



e*5+ 4e ir 5 ir (e i< - 1) 



1 + l + e+ 



9 a+e± , , ?] 



W-hy) 

(3.9) 



up to next-order corrections in the quasi-scale invariant approximation, following what was described 
in Eq. (2.18). These corrections can be absorbed into the coefficients A+, and in table 2. This 
generalizes Eq. (2.20) to rapidly varying interaction vertices. Since (3.9) adds power law contributions 
to the integral, one might worry that not all the previously allowed values of £ will allow convergence of 
the final result - we recall that the overall power law needs to decay faster than (— y) -1 for convergence 
criteria to be met. In the event this does not happen then one is required to perform integration by parts 
an appropriate number of times so as to isolate the primitively divergent contributions, in a completely 
analogous way as we have dealt with the operator £(d Q 2 operator. A convergent answer should similarly 
place constraints on the allowed values of £ such that the correlators do not evolve in time, since this 
could signal a spurious divergence. 

However, for this operator the integral (D.8) has y = 2 — 4e + /(l + e+) and is therefore always 
convergent since the condition Re + £*)] < 3/4 is satisfied for all the range of < e + < 1. The 



16 In DBI models, f x = 1 — cf and therefore one sees that we cannot technically require f x to be constant given that we 
are interested precisely in the regime when the speed of sound is rapidly varying. However, for the purposes of making this 
comparison more transparent, we assume this is the case, similarly to what was done in Refs. [34, 35, 81]. Technically, a more 
precise estimate would have to involve the time dependence of f x in the bispectrum. 
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calculation is carried out in a similar fashion to what we described in 3.1, and which we take as a basis. 
Using the definition of / NL following Eq. (3.2) we find that at leading order in our approximation 



fi 



(leading, equilateral) _ 



NL 



- fc; J (l + eJr(T)cos 



<1 - gj 
1 + g* 



(3.10) 



where T = j— - is a non-negative constant. This result gives the possibility of / NL changing sign if 
e+ > 1/3, when the argument of the trigonometric function changes from the first to the second quad- 
rant. Including the next-order corrections, the result becomes substantially more complicated and in an 
attempt to simplify these expressions as much as possible we evaluate / NL in the so called equilateral 
limit k 1 = k 2 = k 3 = k when k+ = k. 

We organize the correction to the leading order non-gaussianity <5/ NL , in terms proportional to the 
various slow- variation parameters, as follows: 



5/, 
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n s -l 
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(3.11) 
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(3.12c) 



where tn = 2(1 — g*)/(l + g*). We note the dependence on the scale through k, expected whenever the 
slow-roll approximation breaks down. The polygamma functions of order zero, ip^°\ were introduced 
in Eq. (D.4) in appendix D. Also, denotes the m th -harmonic number which relates to i/^ -* via 
^m-i — ip^°\rn) + y E and J satisfies 



j = r(T) 



y E + ln2 + i/> (0) (T) 
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where we have used the results derived in appendix D. Eq. (3.11) is exact at next-order in the quasi 
scale invariant approximation, and each contribution in Eqs. (3.12) to the value of / NL is as important as 
the other. Eqs. (3.10) and (3.11) are to be compared to Eq. (3.15) of Ref. [81]. As pointed out in Ref. 
[34] there is no reason to believe these corrections will be negligible, but their magnitude will depend 
on the values of the parameters e, r\, 5, and n s — 1. We do not attempt to produce an order of magnitude 
of these corrections here, as this is beyond the scope of this paper. We also notice that the dependence 
on the scale k vanishes in the slow-roll limit, when e+ is taken to be a perturbative parameter. This 
observation is in line with our previous comments on the strong scale dependence whenever there was 
an appreciable deviation from the slow-roll regime. 

Burrage et al.. — Finally, as previously discussed, our results extend what was obtained by Burrage et 
al. [34] and also Chen et al. [10] who treated the speed of sound of scalar fluctuations, s, and the 
expansion rate, e, as slowly varying in the P(X,cj>) class of models. Taking the limit of small e and 
restoring s by setting 5 = s + 2e, we find perfect agreement between our results to next-order in slow- 
roll, operator by operator (which is the only sensible way of partitioning the comparison), including the 
logarithmic corrections previously discussed. As emphasized before, these corrections are important to 
correctly evaluate the scale and shape dependence of the bispectrum, as well as in obtaining an estimate 
of non-gaussianities with a certainty level comparable to the resolution of Planck data. As a consequence, 
our results obey Maldacena's factorization theorem [8, 75] in this limit, which represents a non trivial 
consistency check of our calculation. The results presented in Ref. [34] regarding the scale and shape 
dependence of the bispectrum apply to our analysis if the slow-roll approximation is valid. Whenever 
there is significant breaking of the slow-roll approximation, our results diverge. 

4 Conclusions and outlook 

In the present era of precision Cosmology it is crucial to be able to use the experimental data soon to be 
delivered by Planck efficiently. At the same time, the theoretical instruments at our disposal should be 
able to compete with this level of accuracy. The most important of these are the search for non-gaussian 
signatures imprinted in the CMBR which act as a model differentiator. To this end, a primary object of 
study has been the bispectrum of primordial perturbations. 

Whereas most studies of the bispectrum assume strong slow-roll conditions, cosmological con- 
straints still allow for a deviation from this regime, whilst being compatible with inflationary scenarios 
[35] . In this paper we have focused on the phenomenology of inflationary backgrounds where e and s are 
not perturbative parameters, and satisfy the mild requirement that s, \s\ < 1, but combined to produce a 
quasi scale invariant spectrum of scalar perturbations. Under these assumptions we have calculated the 
bispectrum for single-field Horndeski models perturbatively in n s — 1 but to all orders in e and s. We have 
found that both the scale as well as the shape dependence of the bispectrum is encapsulated not only in 
a logarithmic [34], but in a stronger, power-law form. The power-law behaviour is more relevant if the 
breaking from strong slow-roll is stronger. 
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In an optimistic scenario, such behaviour can be used in our attempt to constraint the parameters 
of a given theory more tightly, and potentially eliminate it from the list of sensible models against ob- 
servations. It is quite likely that with Planck the amount of information one is able to extract from the 
CMBR about the early Universe will come to a halt. It is therefore important to be able to retain the scale 
dependence of the bispectrum in our theoretical computations and estimations, which may enable results 
from smaller scales (LSS) [82-85]. In this paper we have focused on obtaining the general expression 
for the bispectrum in Horndeski theories with fixed interaction coefficients, A; . We have also shown how 
to generalize our formulas to include a rapidly varying A, . This work should be regarded as a first step 
towards estimating observables on a background that does not obey the strong slow-roll regime. Because 
these theories have a strong scale dependence this increases our chances to assess the feasibility of these 
theories as working models of the early Universe, by using data from different scales. 

We leave the study of decomposing bispectrum shapes into fundamental harmonics (basis shapes) 
for future work. This partial-wave decomposition is very important in distinguishing between models 
as discussed in Ref. [21]. Also, it will be interesting to discuss the subsequent implications in the 
consistency relations between the amplitudes of each harmonic, which can ultimately allow us to project 
out this class of models based on a dynamics beyond slow-roll. If single field models are excluded by 
observations, we will have to embark on a multi-field inflation exploration, where the interplay between 
curvature and isocurvature modes can be quite intricate (see, for example, Refs. [86-88]) . 
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Table 1. Coefficients of the functions appearing in the B a bispectrum. For simplicity of notation, we define the 
quantities k, Q 1 and Q 2 i n table 7 of appendix E. 
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Table 2. Coefficients of the functions J y appearing in the B a bispectrum, where Cl u and £l 3 + are defined in table 7 
of appendix E. 
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Table 3. Coefficients of the functions appearing in the B h bispectrum for the first three operators, where t?, = 
r-(fc t — 2fc ; ), and S and /(fc l3 fc 2 , fc 3 ) are defined in table 7 in appendix E. 
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Table 4. Coefficients of the functions appearing in the B b bispectrum for the last two operators, where 
£(*t-2k,). 
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operator 



assignments to Jj in B b 



35, 3^71, 3ni 3 
t— — - T - i + fu , 2 -— (n s -l) --(n s -l) 
1 + £* 1 + e* (1 + e*J 4 2 



not applicable 



„ , 4e* 35. Se.rr 3n:i 3 

Cdj^d-H' r- 2 - t— — - n f * + — (n, - 1) -(n s -l 

„ 4e* 35, 3e + Ti + 37ii 3 

d 2 adjd- 2 Cf — - - n f % + — fa - D -(n s -i) 



Table 5. Coefficients of J appearing in the B b bispectrum. 
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Appendices 



A The spectral index beyond exact scale invariance 



In our main formulas we use n s — 1 as a perturbative parameter, but never its explicit formula, since 
this is never necessary throughout the calculation. We have surpassed the need to know this because of 
the special properties of the two-point correlator of single-field models, in particular that it should be 
time independent on super-horizon scales. In this short appendix we present the formula for n s — 1 for 
completeness. In Eq. (2.2) the variable q = a^/zc^ obeys the following differential equation 



dlnq 
dlny 



1 — e —s 



w s 



1 + 



srj + ts 



(1 



From this it follows that 



+ 



2y 2 1 - e - s 2y 2 (1 



1- wx + ts + 



er\ + ts 



sf] 



rp(4- 2s -s +w) 



(A.1) 



where p = 2 + w+s and x = dlnw/diV (which contributes at next-order only in the quasi scale-invariant 
approximation) . We conclude that the spectral index in these theories satisfies 



n s - 1 



3- Wl + 



2p 



+ 



\P 2 



1- wx + ts + 



p(er) + ts) 

1-B -s ' (1-e-s) 2 [_ 2 (1-e-s) 2 
To leading order in the quasi-scale invariant approximation, we find 



(4-2e -s + w) 



(A.2) 



1 = 



— 2e — 3s — w 



(A.3) 



1 — e — s 

in agreement with Ref. [34]. If we instead assume that e and s are strictly constant and consider only 
P{X, <fi) models, we recover the formula deduced in Ref. [35] that 



n c — 1 = — ■ 



2e + s 



(A.4) 



1 — e — s 

which reproduces the exact scale-invariance relation 5 = — 2e. In the slow-roll approximation this reduces 
to 



n s — 1 ~ — 2e — s 



(A.5) 



at leading order in the approximation, which agrees with the results from Ref. [34]. 

In our approximation scheme 5 is a next-order quantity and the only requirement we use in the 
text is that it is defined in such a way that 5 = s + 2s, in P(X,(f>) models. The general formula for 5 in 
all the Horndeski models will depend on the form of w, and therefore x. Using (A.2) we conclude that 



n c — 1 = — 



■|2ex + A(er\ + ts) - etj 



(A.6) 



l + e 3(1 + s) 2 

This formula agrees with our expectation that 5 is indeed a parameter contributing at next-order only. 
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B Dynamics in y - why? 



In this appendix we motivate the use of the time coordinate y instead of the conformal time variable, T. 
Why is this choice of space-time foliation, which defines a hypersurface at each y = constant, preferred 
compared to the more traditional conformal time? Were we to solve for the perturbations starting from 
Eq. (2.1) in t coordinates, we would obtain a formula for the propagator for scalar perturbations, G^, 
which corresponds to the two-point correlator: 

(ax,T)C(y,f)) = G(x,T;y,f) (B.l) 

We proceed to solve the Green's function equation for the propagator, in Fourier space, which reads 
f d 2 Adz 2daW , ,] i 

where k is the comoving wavenumber and the Dirac delta enforces evaluation at t = f. Because scalar 
fluctuations will propagate at a phase velocity c s generically different from that of the light, one usually 
performs a change of variables such that the propagator is a function of the sound horizon z = —kc s r. 
To get the equation for the evolution of in z will demand inverting dz/dr, which without the premiss 
of working with perturbative s becomes algebraically challenging when plugging into the equations of 
motion since s = dc s / dN cannot be treated perturbatively. It turns out that expressing the y evolution of 
background quantities allows to naturally accommodate a rapidly changing speed of sound for the scalar 
perturbations, and therefore large s, avoiding the difficulty just described. In other words, the y variable 
allows to sum all the powers in s. 

Moreover, writing the quadratic action for the fluctuations in the form of Eq. (2.2) makes the 
reproduction of Bessel's equation more transparent, without any need for perturbative expansions. For 
these reasons, our dynamical analysis is presented in y-time, and follows the same lines of analysis of 
Khoury & Piazza in Ref. [35] . 

C Next-order corrections - useful formulas 

When applying the Schwinger-Keldysh formalism we will summon formulas for the corrections of the 
elementary wavefunctions for £ at next order in the quasi-scale invariant approximation scheme. This 
is because one should perform a uniform expansion in the slow-variation parameters which include the 
interaction vertices, but also the propagator for scalar perturbations. In this appendix we collect some of 
the formulas necessary for this expansion. 

Background evolution. — As explained in the main text, we wish to study the evolution of the background 
quantities in the time coordinate y. To do so, we make a Taylor expansion around the time y of horizon 
crossing of some reference scale k* and make a uniform expansion for small 17 * and t*. Such expansion 
is well defined provided we restrict our analysis to a few e-folds after horizon crossing. This ensures that 
e and s have not varied significantly up to that point, such that we can treat 17 and t as perturbative 
parameters in the dynamics. As justified in the main text, this is indeed a good approximation. 
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To get the y -evolution of the speed of sound, c s , we start by evaluating 



dlnc c 



dln(-fc*y) 



1 + 



(i-e*-0 2 



+ 



H-Ky) 



(C.l) 



To avoid cluttering the notation, we introduce the following parameters which contribute at next-order 
only in our approximation 



a = 



P = 



+ 



srj + ts 



1 — e — s ' (1 — e — s) 2 

£7] + tS 

{l-e-sf 



Upon integration, we find 



5. 



/3Jn(-^y)-y(ln(-^y)) 2 



(C.2a) 
(C.2b) 

(C.3) 



1-8+-S+ 

This becomes more simplified if we write = — 2e+ + 5+ and work perturbatively in 5+ (that is, assuming 
0(5*) = 0(tj J = 0(r J), as follows 



c t = c tk (-Kyy+**\i + 



2e. 



1 + 



- 1) 



H-Ky) 



l + e 



{\n{-Ky)Y 



(C.4) 



In the limit when e and s are both strictly constant, thereby resulting in vanishing a* and j8* we reproduce 
the results in Ref. [35] : 



where we have temporarily restored the dependence in s through 5. The additional contributions ap- 
pearing in Eq. (C.4) are precisely the corrections to this purely power-law behaviour and are relevant 
whenever e and s are slowly-varying. 

Proceeding similarly to get the dynamical evolution of the Hubble parameter, we obtain 



H = HX-KyY~ 



i + 



1 - «* - s 

which in the limit of constant e and s reduces to simply 



&ln(-fc*y)-^(ln(-k*y)) 2 



(C.5) 



H ~ (—Ky) 1 -**-** . 

Again, recasting this result in terms of only one non-perturbative parameter, we find 



l + e* 



ln(-fc,y)-y(ln(-fc,y)) : 



(C.6) 



The explicit formulas in Eqs. (C.4) and (C.6) are relevant when replacing in Eq. (2.6) for the scale factor, 
a(y). 

Derivatives of the elementary wavefunctions. — When calculating the bispectra of Horndeski operators 
which are at least once y-differentiated we will require the derivatives of the elementary wavefunctions. 
We find the evolution of the background primordial perturbation is given by 



^background) _ |^*U + g *) 

dy ' k "2^i7(fcc s J 3 / 2 



k 2 ye iky , 



(C.7) 
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whereas the wavefunction corrections to the internal lines obey 

-d< 



d ^internal) _ ^C 1 + g *) , 2 f ..J _ ^ ~ 1 „-ifcy 



+ e ifcy 



5 + e + r?^ 7i n s — 1 

+ i-(n s -l) + ^— ln(-^y) 



1 + (1 + eJ 2 4 



(C.8) 



In obtaining these expressions we have explicitly eliminated the non-perturbative dependence in the 
variation of the speed of sound, parametrized by s, in favour of the perturbative parameter, 5. 

The time variation of the wavefunctions is relevant in correctly reproducing the time independence 
of the correlation functions for £. 

D Useful integrals 

In this appendix we list the integrals which occur when computing the three-point correlators. They 
fall in essentially two different varieties: integrals involving the Ei(£) function (which first appeared in 
Eq. (2.17)) and those which involve power laws and logarithmic integrand functions. We will identify 
master integrals for each of these families and analyse them separately. 

D.l Integrals involving the Ei(^) function 
These integrals appear in the form 

I r (k 3 )= f dy |(-y)r e i(fci+fc2-fc 3 )y f y- e 2ik ^\ (D.l) 

J — OO J —00 ' ' 

where we have explicitly written 7 as a function of the asymmetric momentum (in this case fc 3 ), even 
though it shall be a function of the three momenta, or equivalently k t (perimeter in momentum space, 
k t = ?q + k 2 + k 3 ). The constant y need not be an integer, and this is where the algebra differs from 
that presented in Refs. [17, 34]. In what follows we consider positive, arbitrary y. To simplify this 
integral, one follows the algorithm developed in these references: we apply a transformation of variables 
to convert this integral into its dimensionless version (defining x = i(fc t — 2fc 3 )y) and we make a rotation 
in the complex plane (using w = — ifc 3 ^). Applying the Cauchy integral theorem we arrive at 

^ 3X dw 



w 



e~ 2w } (D.2) 



In this notation 

3 = I \ h and ^3 = r (k t - 2k 3 ) . 
k t - 2k 3 k t 

We will make use of these variables in §3.1. Focusing on the inner integral in (D.2) we note that upon 
integrating by parts and using the expansion series of the exponential function, it has a convergent series 
representation as follows 



rt>3X dw , ^(-l)"(20 3 x) n 

— e- 2w = rE + ln(2g 3 x) + V ' y . (D.3) 



-28- 



7 r (fc 3 ) explicit result 









ln(l-# 3 ) 




#3 + ln(l-# 3 ) 




# 3 (2+# 3 ) + 2 ln(l-# 3 ) 



Table 6. Explicit results for I r integrals with 7 = 0,1,2. As argued in Ref. [34] the function J r (fc 3 ) has no 
singularities in the limit when k 2 , k 3 — » 0, for which # 3 — » and in fact one finds 7 — » — # 3 , Ij — » (—1/2) #3 and 
I 2 — * (—2/3) ^3 for the results above in the table. These precise values guarantee that Maldacena's consistency 
condition [8, 75], which relates the bispectrum evaluated in the squeezed limit and the spectral index, is obeyed. 



Plugging this result into Eq. (D.l) gives 

( 1 Y +1 f r a* 1 ^(-20 3 ) n (n + r )!) 

W = (je-fcj {r(i + r)[r E + in(20 3 ) + ^ o) (i + r)]+l]^ — „„, } , CD.4) 

in which xp^°\z) denotes the polygamma function of order zero and argument z. 17 In obtaining this 
result we used the fact that Re(y) > — 1 for the first two contributions, whereas Re(y) > —2 for the last 
term. This will indeed be true for all the integrals we will analyse for y is strictly non-negative. Finally, 
we can perform the last sum in Eq. (D.4) by observing that it converges for Re(0 3 ) e ] — 1/2; l/2[. We 
conclude that the integral (D.l) has a closed-form representation, given by 

i Y (k 3 )= (^) r+ {r(i + r)[r E + in(20 3 ) + V' (o) (i + r)] 



20 3 r(2 + r) 3F2 (u, 1,2 + r }, {2, 2}, -2G 3 ) I 



(D.5) 



where F is a (convergent) generalized hypergeometric function, HypergeometricPFQ [89] . In the main 
text of this paper we will use the following abbreviated notation 

w= (^k) r+1/ ^ 3) ' (D - 6) 

which, given that I r is a real-valued function, makes the identification of the real part of the result of the 
integral more transparent. Explicit results for positive integers values of y = 0, 1, 2 are listed in table 6. 

D.2 Other integrals 

The other family of integrals which arises in the calculation of the bispectrum is of the form 



dyj^y (-yyU+Bjn{-Ky) + c±(H-Ky)) 2 } , (D.7) 



17 The polygamma function of order m is the (m + l) th derivative of the logarithm of the gamma function, r. We will only 
require polygamma functions of order zero and one in our formulas. 
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where the constant coefficients A*, and C* contain, in general, leading and next-order contributions 
in the quasi-scale invariance approximation. Provided Re(y) > — 1, which is indeed the case for all the 
integrals which will arise in our calculation, this integral converges and gives 



J r = (^) +r r(i + r)R + 



71 



Hkjk t )-i- 



n 



B+ + CJ Hkjk t )-i- 



+ 



71 



b* + 2Q ( HKIK) - i 2 I + c W> (0) (i + r) 



+ V (0) U + r) 
+ cw (1) (i + r) 

Again, to simplify the notation we will refer to this integral in the form 



+ 



(D.8) 



(D.9) 



so as to identify the real contribution to the result. We note that J y is a complex-valued function, but 
given we only require the real part of J Y we need to use Euler's function and write 



1 a i+r 



1 \ 1+r 



cos 



71 

-d+r) 



ism 



71 

-d+r) 



This explains the presence of trigonometric functions in tables 3 and 4. As a final remark, we note that 
these integrals involve logarithmic contributions of the form ln(fc + /fc t ), which are responsible for scale 
dependence, precisely in the same way the power spectrum has a weak logarithmic scale dependence, 
cf. Eq. (2.11). 



E Listing variables 

In the main tables of §3.1, to simplify the formulas, we have introduced a compact notation for combi- 
nations of momenta and slow-variation parameters. This is summarised in the following table. 
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variable 



mathematical expression 



k 


7, i k 2 | k l k 2 k 3 
K t + k t + fc 2 




K 2 


^1^2 + ^1^3 ^2^3 




«1* 


l + 3^+^ln(^) + ^- 


2e,rj, 
U+sJ 2 








^3* 







(n s -l)[-fc t ln2^+2(fc 1 ln2fc 1 /fc t + fc 2 ln2fc 2 /fc t + fc 3 ln2fc3/fc t ) 



/ (fc l5 fc 2) fc 3 ) +3fc(r E - InfcJfcJ - 3^f^ + 4fc t - 2k 2 /k t ] 

n l-i I -I I 4S* I 27IS* 



Table 7. This table collects all the variables used in the tables included in the main text for simplicity of notation. 
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